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Abstract. In this work we study the phenomenological consequences of the existence 
of long-range forces coupled to lepton flavour numbers in solar neutrino oscillations. 
We study electronic forces mediated by scalar, vector or tensor neutral bosons and 
analyze their effect on the propagation of solar neutrinos as a function of the force 
strength and range. Under the assumption of one mass scale dominance, we perform 
a global analysis of solar and KamLAND neutrino data which depends on the two 
standard oscillation parameters, Am^i and tan^0i2, the force coupling constant, its 
range and, for the case of scalar-mediated interactions, on the neutrino mass scale as 
well. We find that, generically, the inclusion of the new interaction does not lead to 
a very statistically significant improvement on the description of the data in the most 
favored MSW LMA (or LMA-I) region. It does, however, substantially improve the fit 
in the high-Am^ LMA (or LMA-II) region which can be allowed for vector and scalar 
lepto-forces (in this last case if neutrinos are very hierarchical) at 2.5 a. Conversely, 
the analysis allows us to place stringent constraints on the strength versus range of the 
leptonic interaction. 
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Many extensions of the standard model of particle physics predict the existence of 
new forces; either generating deviations of the gravitational law at short distances or 
predicting low mass particles whose exchange will induce new forces at long distances, 
generally violating the equivalence principle I2]- A number of experiments have 
been searching for such new forces but no evidence has yet been found. These null 
results, nevertheless, provide very significant constraints on particle physics models, 
gravitational physics and even on cosmological speculations 

We know that there are two long-range forces in nature, namely the electromagnetic 
and the gravitational force, but there is no a priori reason why only these two should 
exist. Since the seminal work of Lee and Yang |1] it has become standard to consider 
long-range forces coupling to baryon and/or lepton number. These would definitely 
violate the universality of free fall and thus can be tested by Eotvos-type experiments, 
as pointed out in [4]. In particular, Okun used this idea to establish a bound on 
the strength of an hypothetical vectorial leptonic force; see ^ for a recent review. He 
obtained the bound on the "fine structure" constant 

kv < 10"^9 . (1) 

The lepton flavour symmetries, as several solar, atmospheric, as well as KamLAND, 
K2K and MINOS neutrino experiments indicate, cannot be exact in nature. If an 
electronic force exists, we may expect it to be of arbitrary but finite range. Notice that 
when the range is less than the Earth-Sun distance, the bound (0) is no longer valid. 
Other experiments |2l E] using the Earth instead of the Sun as the electronic source 
derive bounds, which, however, are much less strict than (^. 

Recently, it has been shown that neutrino physics allow to place quite stringent 
bounds on these new interactions. The impact of vector-like long-range forces coupled 
to lepton flavour numbers on solar and atmospheric neutrino experiments has been 
investigated in [7j and [S], respectively. It has been estimated in p| that the fine structure 
constant coupled to electronic number ky{e) < 6.4 x 10~^^ to be compatible with 
solar neutrino oscillations. According to 8j atmospheric neutrino data can constrain 
the strength of forces coupled to — L^, kv{efi) < 5.5 x 10~^^, and to Lg — L^., 
kvier) < 6.4 x 10^^^ (at 90% CL) when the range of the force is the Earth-Sun 
distance. Also, new forces coupling individually to the muonic or tauonic number 
can be constrained, producing, however, weaker bounds. In particular, primordial 
nucleosynthesis considerations provide the bound A;y(/i, r) < 1.8 x 10^^^ P]. 
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In the present paper, we would like to continue the line of research started in 0. We 
perform a global analysis of the solar and KamLAND neutrino data with the objective 
of constraining new electronic forces. Since long-range forces can be mediated by the 
exchange of massless or very light bosons of spin J = 0, 1 and 2, we study all these 
cases. In the next Section the formalism to take into account the effects of scalar (S), 
vector (V) and tensor (T) forces in neutrino oscillations is worked out. In Section El we 
perform the analysis and we devote the final Section to our conclusions. 

2. Formalism: Effects in Solar Neutrino Oscillations 

If there is a new force coupled to the leptonic flavour numbers, its presence will affect 
neutrino oscillations when neutrinos travel through regions where a flavour dependent 
density of leptons is present as it is the case for solar neutrinos. 

As we will see more in detail, how the effect of the new interaction on the oscillation 
pattern depends on its Lorentz structure. Nevertheless, it can be easily proved that 
for scalar, vector or tensor interactions of large enough range the modification of the 
evolution equation can always be casted in terms of a unique function which solely 
depends on the background density of leptons - the source of the force - and the range 
of the interaction. 

In particular, in the Sun, we can define a universal function W determining the 
effect of a new force coupled uniquely to Le at a point r from the center of the Sun: 



where -Rsun is the radius of the Sun and ne(r) is the electron number density in the 
medium. The range of the interaction is A = 1/m, where m is the exchanged boson 
mass. We can easily do the angular integration in Q and get 



In Figure n we show the function W{r) in the Sun as a function of the distance in 
units of -Rsun for various ranges A. As seen in the figure, inside the Sun and up to 0.1 
-Rsun, W is constant independently of the range of the force considered. After that, it 
drops according to the force range and at the position of the Earth, at about 200 -Rsun 
from the center of the Sun, W < 2 x 10^^^ eV. 

In order to see the effect of the new Lg-coupled force in the evolution of solar 
neutrinos we have to take into account that at present the minimum joint description of 
atmospheric K2K MINOS solar [13 IH HSl CEl HIl EE! and reactor [IHlllI! 




(2) 
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data requires that all the three known neutrinos take part in the oscillations. The mixing 
parameters are encoded in the 3x3 lepton mixing matrix which can be conveniently 
parametrized in the standard form 
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where Cj^ = cos 9ij and s^- = sin^^j^. 

According to the current data, the neutrino mass squared differences can be chosen 
so that 



|Am^2l 



Am? 



(5) 



As a consequence of the fact that Am2i/\Amli\ ^ 0.03, for solar and KamLAND 
neutrinos, the oscillations with the atmospheric oscillation length are completely 
averaged and the interpretation of these data in the neutrino oscillation framework 
depends mostly on Amgi, 612 and 613 (while atmospheric and long baseline neutrino 
oscillations are controlled by Am|2, 6*23 and ^13). Furthermore, the negative results 
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from the CHOOZ reactor experiment 20 imply that the mixing angle connecting the 
solar and atmospheric oscillation channels, 6'i3, is severely constrained (sin^ ^13 < 0.041 
at 3(T Altogether, it is found that the 3-z/ oscillations effectively factorize into 2-z/ 

oscillations of the two different subsystems: solar (and reactor) and atmospheric (and 
long baseline). 

Because the new interaction is flavour diagonal and its effect in the evolution of 
atmospheric neutrinos does not modify the hierarchy (0), the 21/ oscillation factorization 
still holds. Thus, in general, one can write the solar neutrino evolution equation in the 
presence of the new Lg-coupled force as 

' V{r) 








(6) 



where 



/ cos 6*21 sin 6112 \ 
Y -sme^i2 cos6'i2 J 

and M and V{r) will depend on the Lorentz structure of the new leptonic interaction 
as we describe next: 

2.1. Scalar Interaction (J = 0) 

If the new force is mediated by a neutral spin J = particle, 0, the interaction 
Lagrangian for the neutrinos will take the form: 

L = -go(t)i'uipu ■ (8) 

We now use that the electron density can be taken to be the (static) source of the Lg- 
coupled interaction and that the neutrino will experience the force due to this source as 
long as within its range. In this approximation the effect of the new interaction can be 
included in the neutrino evolution equation as in (jHI) with: 

M= ' n n Vir) = Vccir), (9) 

where mi 2 are the neutrino masses in vacuum, Vcc{i^) = V2GFne{r) is the Mikheyev- 
Smirnov-Wolfenstein (MSW) potential and Ms(r) = ks{e)W{r) is the new scalar- 

leptonic- force-induced mass term. Here ks{e) = — ^. Different from the vector and 
tensor cases, the z/g survival probability will not only depend on Arri^i = rn^ — m\ but 
also on the value of the absolute neutrino mass scale mi, which is still not experimentally 
known [22] ■ Notice also that the term Ms{r) has the same sign for neutrinos and 
antineutrinos. 
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2.2. Vector Interaction (J = 1) 

If there is a new Lg-coupled vector force mediated by a neutral vector boson Aa, the 
interaction Lagrangian for the neutrinos is: 

L = -giAj.-f^^^. . (10) 

In this case the neutrino evolution is described by Q with 

M=('"' M and V{r) = Vcc{r) + VY{r), (11) 
\ m2 J 

where Vy = kv{e)W{r) is the lepto- vector potential with /cy(e) = As a consequence 

Air 

of the vector structure of the force, the new leptonic potential adds to the MSW potential 
with the same energy dependence and sign. It will accordingly flip sign when describing 
antineutrino oscillations. 

2.3. Tensor Interaction (J = 2) 

For a Le-coupled force mediated by a tensor field of spin J = 2, Xap, the interaction 
Lagrangian for the neutrinos is given by 

L = -^ Xap [ipur^d^tP, - za"^,7^V'.] , (12) 

where the coupling constant g2 has dimensions 1/E. 

In this case the neutrino evolution is described by (jB)) with 

M = I ^ ] and V(r) = Vccir) - VT(r). (13) 

V m2 y 

Vrlr) = Ei,kT{e)W{r) is the tensorial potential and is the neutrino energy. Here 

krie) = rUe— and rUe is the mass of the electron. 

Finally, let us remark that, as for the case of a scalar leptonic force, a tensor force is 
always symmetric when changing from neutrinos to antineutrinos. This is obvious from 
fll2|) since what we have coupled to the tensor field x is in fact the energy momentum 
tensor of the leptons, which has to be symmetric under the exchange of particles and 
antiparticles. 

Before performing the data analysis we can have some idea of the size of the effects 
of these long-range forces on neutrino oscillations by studying the value of the survival 
probability of solar neutrinos for different values of the coupling constant and range. The 
probability is obtained by solving © along the neutrino trajectory from its production 
point in the Sun to its detection point on the Earth. We have verified that for the 
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Figure 2. Survival probability of in the Sun as a function of the neutrino energy 
for an infinite range scalar (first two panels), vector (third panel) and tensor (lower 
panel) ig-coupled force, for various values of the strength fci(e), i ~ S,V,T. For all 
curves we have used tan^ 9i2 — 0.44 and Am^i — 7.9 x 10^^ eV^. 



range of parameters of interest the evolution in the Sun and from the Sun to the Earth 
is always adiabatic. It is interesting to remark that for the largest range forces the 
effective mixing angle at the surface of the Earth is still affected by the value of the 
leptonic potential produced by the Sun electron density. 

We show in Figure |21 the survival probability of solar z/g at the sunny face of the 
Earth as a function of the neutrino energy E,^ for an infinite range force. In the upper 
panels we show the scalar case for two extreme values of mi, which correspond to strong 
mass hierarchy (mi = 0) and degenerate spectrum (mi = 0.1 eV), and some values of 
ks{e). From this figure we can foresee that for mi = 0, values of ks{e) > 10~^^-10~^^ 
will conflict with the existing solar neutrino data while for mi = 0.1 eV even smaller 
values of the coupling, ^^(e) > 10~^^-10~^^, will be ruled out. In the third panel we 
show the vector case for some values of kv{e). One expects from this that our analysis 
will lead to bounds kv{e) < 10~^'^-10~^^. Finally, in the lower panel the tensor case is 
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Figure 3. Same as FigureHlbut for A/i?sun = 0.01,0.1, and 1. 



displayed. In this case one expects that the data will constrain kx[e) < 10 — 10 
eV-^ 

The dependence of the v^. survival probability on the range A/-Rsun can seen in 
Figure 01 From this plot we can see how the limit on k will weaken as we allow for the 
force to have a finite range. 

3. Constraints from Solar and Reactor Neutrino Data 

We present in this section the results of the global analysis of solar and KamLAND data 
for the long-range forces discussed in the previous section. 

Details of our solar neutrino analysis have been described in previous papers [221 1211 • 
We use the solar fluxes from Bahcall et al [23 • The solar neutrino data includes a total 
of 119 data points: the Gallium P^ IT3] (2 data points) and Chlorine jT3| (1 data point) 
radiochemical rates, the Super-Kamiokande JHl zenith spectrum (44 bins), and SNO 
data reported for phase 1 and phase 2. The SNO data used consists of the total day- 
night spectrum measured in the pure D2O (SNO-I) phase (34 data points) [T71, plus 
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the full data set corresponding to the Salt Phase (SNO-II) jTHj. This last one includes 
the NC and ES event rates during the day and during the night (4 data points), and 
the CC day-night spectral data (34 data points). It is done by a analysis using the 
experimental systematic and statistical uncertainties and their correlations presented 
in JH], together with the theoretical uncertainties. In combining with the SNO-I data, 
only the theoretical uncertainties are assumed to be correlated between the two phases. 
The experimental systematics errors are considered to be uncorrelated between both 
phases. 

In the analysis of KamLAND, we neglect the effect of the long range forces due 
to the small Earth-crust electron density in the evolution of the reactor antineutrinos 
and we directly adapt the map as given by the KamLAND collaboration for their 
unbinned rate+shape analysis j_26^ which uses 258 observed neutrino candidate events 
and gives, for the standard oscillation analysis, a Xmin~'^^^-'^^- The corresponding 
Baker-Cousins for the 13 energy bin analysis is Xmm — 13.1/11 dof. 

3.1. Scalar Interaction 

In the scalar case the analysis of solar and KamLAND data depends on 6 parameters: 
the two oscillation parameters, Arri^i and tan^ 612, the absolute neutrino mass scale mi, 
the coupling ks{e) and the range of the force A/-Rsun- We have done our analysis fixing 
nil at two extreme values, mi = and 0.1 eV. 

For mi = eV (0.1 eV) we find that the best fit point is 

tan^ 612 = 0.49 (0.44) Am^^ = 7.9 x 10"^ eV^ (7.9 x 10^^ eV^) 
\/Rsun = 1.26 (indet) ksie) = 5 x 10"^^ (0) 

AxLn = -1-4 (0) , (14) 

where the range of the interaction is indeterminate (indet) when ks{e) = and A^^j^ 
is the shift in compared to the best fit point in the absence of a new Lg-couple force 
for which 

tan^ 9i2 = 0.44 and Am^^ = 7.9 x 10"^ eV^ . (15) 

Thus we find that for highly hierarchical neutrinos the presence of the long range force 
can lead to a small improvement of the quality of the fit but this improvement is not 
statistically very significant leading only to a decrease of less than two units in even at 
the cost of introducing two new parameters. For degenerate neutrinos, on the contrary, 
the presence of the long-range force does not improve the agreement with the data. 
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Figure 4. Allowed regions from the global analysis of solar plus KamLAND data in 
the (Am2i, tan^ ^12) parameter space at 3cr CL (2dof) for the different type of long- 
range forces as labeled in the figure. The standard MSW allowed region (shaded area) 
is also shown for reference. 



This is further illustrated in Figure E] where we show the result of the global analysis 
of solar data plus KamLAND data in the form of the allowed two-dimensional regions 
at 3cr CL in the (Am2i, tan^ ^12) plane after marginalization over ks{e) and A. The 
standard MSW allowed region is also showed for reference. As seen in the figure for 
any value of mi the presence of the long-range forces has a mild (or null) impact on the 
allowed range of Am2i and tan^ 612 in the low-Am^ island, known as large mixing angle 
region I (LMA-I). Most interestingly, we also find that in the presence of the long-range 
scalar force the description of the solar data in the high-Am^ (LMA-II) region can be 
significantly improved if neutrinos are very hierarchical. As a consequence, in this case, 
there is a new allowed solution at the 98.9% CL. The best fit point in this region is 
obtained for 

tan^ 9i2 = 0.5 Am^^ = 1.6 x 10"^ eV^ 

A/i?sun = 4.0 ksie) = 6.9 x 10"^^ 

AxLn = 7.2 (16) 
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While this region is excluded at more than 4ct for the standard MSW oscillations it is 
allowed at 2.5cr in the presence of a new scalar force for a narrow band of the lepto- 
force parameters. For example, for A = oo the LMA-II region is allowed as long as 
4.5 < {ks{e)/10~^^) < 8.0 t 

The improvement of the fit in the LMA-II region can be understood as follows. For 
the range of scalar Le-coupled lepto-forces that we are considering the neutrino evolution 
in the Sun is always adiabatic. Therefore the corresponding survival probability is 
determined by the effective mixing angle at the neutrino production point, tq, which 
takes the form 

S2i2^m = — I ' (17) 



'[Aml sC2u - 2E,Vcc{ro) - M|(ro) + Ms{ro){rm + m^W + S2^2Aml^s 

where we have defined Ani^i g = Arri^i — Msir^Am^x with Am^x = m2 — mi and 
5*212 = sin26'i2, and C2i2 = cos2^^i2. As long as the term Ms{rQ){mi + m2) in the 
denominator is subdominant, the net effect of MsItq) is to increase the effective matter 
potential and correspondingly, one can achieve a good fit to the solar data with higher 
values of Arri^i. In this case, the CL at which the LMA-II region is allowed is determined 
by KamLAND data ^19^ because the fit to the solar data cannot discriminate between the 
LMA-I and LMA-II regions once scalar Lg-coupled force effects are included. Clearly 
this implies that this solution will be further tested in the future by a more precise 
determination of the antineutrino spectrum at KamLAND. 

As mi increases the term Ms{ro){mi + 7712) becomes more important and 
correspondingly the quality of the fit within the LMA-II region worsens. We find that 
for nil ^ 0.02 eV there is no LMA-II solution at 3a. 

We show in Figure El the survival probability for this best fit point in the high-Am^ 
(LMA-II) region in the presence of the scalar lepto-force together with the extracted 
average survival probabilities for the low energy {pp) , intermediate energy (^Be, pep 
and CNO) and high energy solar neutrinos (^B and hep) from [2H1- For comparison we 
also show the survival probability for conventional MSW oscillations for the same values 
of Am^i and ^12. From the figure it is clear that the inclusion of the lepto-force, leads 
to an improvement on the description of the solar data for all the energies being this 
more significant for intermediate- and high-energy neutrinos. 

Conversely, the global analysis of solar and KamLAND data results into constraints 
on the possible size of the new contribution to the neutrino mass. This is illustrated 
in Figure IHl where we show the result of the global analysis in the form of the allowed 
I The LMA-II region was found also for generic environmentally dependent solar neutrino masses |27| 
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Figure 5. survival probability in the Sun versus neutrino energy for the best fit 
point in the high-Am^ region in the presence of scalar lepto-force with parameters 
given in (|16|l (full line) and for the vector lepto-force with parameters in 121|l . The 
dotted line is the survival probability for conventional MSW oscillations {h = 0) for 
the same values of Am^x and 9i2. The data points are the extracted average survival 
probabilities for the low energy (pp), intermediate energy (^Be, pep and CNO) and 
high energy solar neutrinos (^B and hep) from j28| . 



two-dimensional regions in the (^^(e), A/-Rsun) parameter space after marginalization 
over Am^i, tan^ 612 for the two reference values of mi. 

As seen in the figure, already at 90% CL (or lower) the regions are connected to 
the standard ks{e) =0, X = 00 case. In other words, the analysis shows no evidence 
of any long-range Lg-coupled scalar force contributing to the solar neutrino evolution. 
However, as expected, there is a strong correlation between the implied bound on the 
strength and the range of the interaction. As discussed in the previous section the 
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Figure 6. Allowed regions from the global analysis of solar plus KamLAND data in 
the fcs(e) versus A/i?sun plane at 90%, 95%, 99% and Sa CL (2dof) for a ie-coupled 
scalar force, with mi = eV (left) and mi = 0.1 eV (right). 



tightest bound on the strength of the new interaction corresponds to A = oo (which 
effectively apphes for A > 10 -Rsun) 

ks{e) < 5.0 X 10~^^ (^0 < 2.5 x lO^^^) for mi = eV , (18) 
ks{e) < 1.5 X 10~^^ (^0 < 4.3 x lO^^a^ ^ ^ ^^g^ 

at 3a (Idof). 

3.2. Vector Interaction 

If the new leptonic interaction is mediated by a vector boson the analysis of solar and 
KamLAND data depends on 4 parameters: the two oscillation parameters, Am^i and 
tan^6'i2, the coupling kv{e) and the range of the force X/Rsnn- 
In this case the best fit point corresponds to 

tan^ 012 = 0.46 Am^i = 7.9 x 10"^ eV^ 

\/Rsnn = 2.5 X 10"^ kv{e) = 7.9 x 10"^^ 

AxLn = -0.8 , (20) 

which as for the scalar case represent a very small variation with respect to the standard 
LMA solution as seen in Figure lU 

From Figure |3 we also find that the presence of the vector lepto-force can improve 
the quality of the fit in the LMA-II region. This is possible because the lepto-force also 
adds up to the MSW potential and consequently the vacuum-matter transition condition 
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Am2i cos 2^ = 2E^V{rQ), which has to occur around Ey 1 MeV to fit the data, can 
be achieved for a larger value of Am^i. 

So also in the vector case, there is a new allowed solution at 2.5 a CL. The best fit 
point in this region is obtained for 

tan^ 612 = 0.46 Amj^ = 1.6 x 10"^ eV^ 

A/i?sun = 4.0 kvie) = 6.9 x 10"'*^ 

AxLn = 7.9 (21) 

The corresponding survival probability for this point is shown in Figure El The LMA-II 
region is allowed for a narrow band of {kv{e), X/Rsnn) values. For example for A = 00 
the LMA-II region is allowed as long as 2.0 < (A;v(e)/10~^^) < 30. 

As for the scalar lepto-force with hierarchical neutrinos, the CL at which the LMA- 
II region is allowed is determined by KamLAND data [19^ and consequently this solution 
will be further tested by this experiment. 

Concerning the allowed ranges of the lepto-force parameters, we show in Figure [7| 
the result of the global analysis in the form of the allowed two-dimensional regions in the 
{kv{e), A/i?sun) parameter space after marginalization over Amgi, tan^ 612. In particular 
we find that for A = cxd 

kv{e) < 2.5 X 10"^^ (^1 < 1.8 x 10"^^) , (22) 

at 3(7 (Idof). 

3.3. J = 2 Interaction 

In the tensor case the analysis of solar and KamLAND data depends on 4 parameters: 
the two oscillation parameters, Aml^ and tan^6'i2, the coupling /cr(e) and the range of 
the force A/-Rsun- 

We find the best fit point at 

tan^ 9i2 = 0.44 Am^^ = 7.9 x 10"^ eV^ 

A/i?su„ = (indet) krie) = 

AxLn = . (23) 

Thus for the case of the J = 2 interaction there is no improvement on the quality of the 
fit compared to the standard MSW scenario. 

From Figure |3] we see that the only effect of the inclusion of this interaction is the 
enlargement of the LMA-I region for smaller values of the mixing angle. The inclusion 
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Sun 

Figure 7. Allowed regions from the global analysis of solar plus KamLAND data in 
the fcy(e) versus A/i?sun plane at 90%, 95%, 99% and 3a CL (2dof) for a Le-coupled 
vector force. 

of the J = 2 lepto- force decreases the effective matter potential (see (IT!^ ) at higher 
energies. As a consequence, as seen in Figures |2] and El the corresponding survival 
probability grows for the high end of the solar neutrino spectrum as compared to the 
standard MSW case. This behaviour can be compensated by a decrease of the vacuum 
mixing angle and leads to the enlargement of the allowed LMA-I region. However, 
because the effect of the J = 2 lepto-force is to decrease the matter potential, no LMA- 
II solution is found in this case. 

Finally, in Figure |Hl we show the result of the global analysis in the form 
of the allowed two-dimensional regions in the (^^(e), A//2sun) parameter space after 
marginalization over Amgi and tan^ 612. For A = 00 we obtain the bound 

Me) < 1.7 X 10-^° eV~^ (^2 < 6.5 x 10"^=^ eV"^) , (24) 

at 3a (Idof). 

4. Conclusion and discussion 

In this paper be have investigated the effects of long-range Lg-coupled forces in the 
oscillation of neutrinos in the Sun. Their impact on neutrinos depend on the spin of the 
exchanged particle leading to the new force. We have treated, in turn, the consequences 
of scalar, vector and tensor forces. 
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?i/R„ 

Sun 

Figure 8. Allowed regions from the global analysis of solar plus KamLAND data in 
the kxi^e) versus A/i?sun plane at 90%, 95%, 99% and 3a CL (2dof) for a ie-coupled 
tensor force. 

In our study we have used data from solar neutrino experiments and KamLAND. 
In the global fit, apart from the mass difference Am2i and the mixing angle 612, which 
are the usual standard parameters, we have the new physics parameters ki, the strength 
of the interaction with i=S,V,T, and the range A of the interaction. For the scalar case, 
we also need to assume the absolute mass scale for neutrinos. 

The result of the global fit is that there is not a significant improvement in the 
description of the data in the preferred (in the standard case) MSW LMA region (LMA- 
I). However, the addition of new physics permits to "recover" the high-Am^ region 
(LMA-II) which is now allowed in the case of scalar and vector forces. Why this is so 
only for these forces and not for the tensorial one can be understood quite simply, as 
we have explained in the paper. 

We have deduced limits on the strength of the forces for different values of the 
range A (see Figures IHl 13 and IB)) . For an infinite range (A = 00) Lg-coupled force, we 
get the following limits in terms of the "fine structure constants" ki 

ks{e) < 5.0 X 10"^^ (mi = eV) , (25) 
ks{e) < 1.5 X 10"^^ (mi = 0.1 eV) , (26) 
fcv/(e) < 2.5 X 10-^^ (27) 
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krie) <1.7 xlQ-^^ eV-\ (28) 

at 3cr (Idof). 

These bounds are practically the same for any A > 10-Rsun- For A < 10-Rsun, the 
bound slowly worsens, and for A ~ 0.1i?sun we have that the bound on ki is a factor 
less than 10 weaker than ()25|) - (f^ . In |7j the vectorial case was evaluated, not with the 
statistical rigor we have proceeded in the present paper. A bound on kv{e) one order 
of magnitude more stringent was found. Apart from extending the results to include 
scalar and tensor forces, our bound ^I7\ supersedes the one in [?]. 

In deriving ()25j) to ()28j) we have assumed that the new force only couples to Le 
but the analysis can be easily extended to forces coupled to any linear combination 
of the lepton flavour numbers Lg + aL^ + For vector and tensor forces it can 

be easily proved that as long as the 2z/-oscillation factorization holds, the bounds 
will be the ones given in ()27|1 and but they will constrain the combination 

kv,T = A;v'r(e)(l — a cos 2^23 sin 26'23). In particular this means that for the anomaly- 
free combinations L^, — L^, and Lg — Lj. the bounds on the corresponding ky^J■{^^J') and 
kv,T{GT) become tighter by a factor 1 + cos2^^23 — 2 and 1 + sin2^^23 — 2 respectively. 
For the case of the scalar interactions coupled to a general linear combination of lepton 
flavour numbers, the combination of couplings bounded depends on mi and it does 
not have such a simple expression as for the tensor and vector cases. But generically, 
provided that there is not a fine-tuned cancellation between the masses and the new 
interaction terms, the bound will be of the same order of magnitude as given in ()25p 
and (j^. 

A final issue we should care about is the possibility of screening of leptonic charges 
by electronic neutrinos and antineutrinos from the cosmological relict background. As 
is explicitly demonstrated in [7], screening plays no role at the relevant scales on the 
order of the solar radius. 
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